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NOMENCLATURE 



rotor blacle chord 

rotor thrust coefficient 

Rotor lag damping coefficient 

rotor blade flapping moment of Inertia 

rotor blade pitch moment of Inertia 

rotor radius 

time to half amplitude of mode (seconds) 
time to double amplitude of mode (seconds) 
period of mode (seconds) 
forward speed 

distance blade section aerodynamic center Is aft of the 
pitch axis (fraction of chord) 

distance blade section center of gravity Is aft of the 
pitch axis (fraction of chord) 


/- 


to® 

XI 


blade Lock number 
damping ratio of a root 

eigenvalue or root of a mode (dimensionless, based on the 
rotor rotational speed) 

zocor advance ratio, forward speed divided by tip speed 
flap natural frequency In rotating frame (per rev) 
blade lag TOtatlng natuzal frequency (per rev) 
rotor solidity ratio 

natural frequency of blade pitch motion (per rev) 
rotor rotational speed 


ill 


ELEMENTARY APPLICATIONS OF A ROTORCRAFT 


DYNAMIC STABILITY ANALYSIS 

Wayne Johnson* 

Ames Research Center 
and 

Ames Directorate 

U.S. Army Air Mobility T &D Laboratory 


SUMMARY 

A number of applications of a rotorcraft aeroelastic analysis are 
presented, intended to verify that the analysis encompasses the classical 
soliit1nn« of rotor dynamics, and to examine the influence of certain 
features of the model, Results are given for the following topics j flapping 
frequency response to pitch control? forward flight flapping stability; 
pitch/flap flutter and divergence? ground resonance instability? and the 
flight dynamics of several representative helicopters. 

INTRODUCTION 

An aeroelastic analysis has been developed for rotorcraft, applicable 
to isolated rotor, rotor with wind tunnel support, and complete helicopter 
dynamics. The analysis is described in reference 1, This report presents 
a number of examples of the dynamic behavior calculated using this aeroelastic 
model. The first purpose of these applications is to check out the 
sophisticated analysis by demonstrating that it encompasses several elementary, 
classical solutions in rotor dynamics. The second purpose is to examine 
the influence on the dynamics of some features available in the model, 

RESULTS AND DISCUSSION 

Flapping Frequency Response 

Figures 1 to 3 present the frequency response of a threq-bladed rotor 
in hover to pitch control inputs. Figure 1 shows the coning response to 
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collective inputn for an articulated rotor ( ^ •« 1 and P). Figures ?. and 
3 show respectively the lateral and longitudinal tip-path- plane tilt response 
to lateml cyclic input for a hingeless rotor "' 1.15 and * B), Three 
cases are shown*, with no inflow lynatnics* with a quaBlstatlo inflow dynamics 
model, and with the complete inflow dynamics model. The difference between the 
quasistatic and complete inflow dynamics models is simply that the latter includes 
a first order time lag in the inflow response to rotor force perturbations (see 
reference l), The influence of the rotor Inflow on both the steady and high 
frequency resi»nse can be significant, but the quasistatic model gives nearly the 
sane results as the complete inflow dynamics model (c.f. references 2 and 3). 

Forward Flight Flapping Stability 

Figure presents a root locus of the flapping stability of an 
articulated rotor blade in forward flight « 1 and » lO). The blade 
motion consists of the rigid flap degree of freedom, with no shaft motion. The 
analysis is for a single blade in the rotating system, including the influence 
of the periodic coefficients on the dynamic stability. This result, including 
the advance ratio at the stability boundary, compares well with simpler solutions 
(c.f, references and 5 )» The roots were calculated using a step size of 

12° in the integration of the equations of mo tion (over one period , for 
use in the Floquet theory analysis) , This azimuth increment is really more 
appropriate for a three- or four-bladed rotor; for a single blade at high 
af^vance ratio a step size of -oH’ =* 1 or 2° gives more accurate results . 

Pitch/F'lap Flutter 

Next consider the pitch/flap flutter of an articulated rotor blade 
(<7 =1 and ^-8). The degrees of freedom involved are rigid flap and 
rigid pitch motion, with the pitch bearing outboard of the flap hinge. 

The following va,lues are used for the rotor chord (fraction of rotor radius), 
pitch inertia (fraction of flapping Inertia), and rotor thrust coefficient 
to solidity ratio: c/H “ 0.1, “ 0,001, and 0.^®^ “ 0.0?. The 

flutter and divergence boundaries are plotted in figures 5 to 7, as the 
control system stiffness (l.e. the pitch natural frequency lOq ) required 
for stability as a function of the chordwise center of gravity/aerodynamic 
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center offsot# A large enough control system stiffness or forward center 
of gravity position stabilises both the divergence and flutter modes# 

Figure 5 preaentn the pltch/flap flutter boundaries for the hovering 
articulated rotor with aerodynamic center offsets " 0 and ,05o# This 
result compares well with the classical solutions (see reference 6). 

Figure 6 shows the affect of adding the first flapwise bending mode 
(with natural frequency ?.6/rev) to the analysis, for the hovering rotor 
vfith • 0. Finally figure ? shows the influence of forward flight on 
the flutter boundaries, for advance ratios /a •* 0 and ^ « 0.3* The 

.3 results are based on a constant coefficient approximation, using 
the average of the coefficients in the rotating frame} it is better to make 
this approximation in the nonrotating frame however, for a rotor with 
several blades# 

Ground Resonance 

Now we shall consider some results from ground resonance investigations. 
Ground resonance is a mechanical instability involving the inplane motion 
of the rotor hub due to flexibility of the helicopter or wind tunnel support, 
coupled with the rotor cyclic lag degrees of freedom (see reference ?)# 

For a specific i’otor support system, the ground resonance stability boundary 
may be expressed in terms of the lag damping required to stabilize the motion, 
Beutsch (reference 8) developed a simple criterion for the lag damping 
required for stability, based on the assumption that the blade mass is 
small compared with the support mass (which is almost always true). 

Figure 8 presents a comparison of the ground resonance stability boundary 
calculated using the dynamics analysis of reference 1, with the results of 
the Deutsch criterion. The points are for two full-scale rotors on a 
number of wind tunnel support configurations; the line represents exact 
agreement of the two calculations. The correlation is very good. 

Figure 9 presents the ground resonance stability calculated for a 
hovering, four-bladed articulated rotor in a wind tunnel. The rotor and 
the support properties for this case are given in table 1 . The analysis 
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used the rigid flap and lag degrees of freedom for each blade » and two 
longitudinal and two lateral degrees of freedom for the hub support. The 
figure shows the damping ratio of the four support modes as a function 
of rotor speed , for three levels of blade lag damping (G ^ ) . Each support 
mode is stabilized at the resonance with the high freq.uency lag mode — 
^support ** )S1{ arid at a higher rotor speed a degradation in stability 

occurs due to the renonarice with the low freq.uency lag mode — *^support *" 

(l - With low enough lag damping it is the latter resonance which 

produces the ground resonance iristability. In general the mode which goes 
unstable in this resoriance (l.e. the support mode or the lag mode) depends 
on which has the least uncoupled damping. In the root locus it is the 
mode nearest the imaginary axis which actually makes the excursion into 
the right-half-plane as the rotor speed is varied. In the example of figure 
9 , the lateral balance mode has a low level of damping (table l), so at 
resonance the support mode damping decreases while there is a corresponding 
increase in the damping of the low frequency lag mode. For well damped 
support modes, the ground resonance instability will occur at such low 
lag damping that it is the lag mode which actually goes into the right- 
half-plane (a mode involving both rotor and support motion of course). 


Helicopter Flight Dynamics 


The rotorcraft aeroelastlc analysis of reference 1 has also been 
applied to the calculation of helicopter flight dynamics. Tables 1 to 5 
present typical results. Four aircraft are considered! a small articulated 
rotor helicopter (gross weight Il60 kg, rotor radius 4.0 m), a large 
articulated rotor helicopter (gross weight 15200 kg, rotor radiiAs 11.0 m), 
a soft-inplane hingeless rotor helicopter (gross weight 2100 kg, rotor 
radlv© 4.9 m), and a tandem rotor helicopter (gross weight B600 kg, rotor 
radius ?.6 m). The tables give the period and damping for the flight 
dynamics modes in hover and in forward flight at 100 knots. The basic 
model has three rigid body degrees of freedom (uncoupled lateral or longitudinal 
motions) plus the rotor flapping degrees of freedom. The roots for the 
basic model are compared with results using a quasistatic approximation 
for the rotor flap dynamics? including the rotor lag degrees of freedom; 
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Including the rotor inflow perturhaiion (unsteady aerodynamics)! and with 
results including all six rigid body degrees of freedom (coupled longitudinal 
and lateral dynamics )• For the soft^inplane hingeless rotor (table 4), 
the results Including both rotor lag and rotor torsion degrees of freedom 
are also presented. For the tandem rotor helicopter (table 5)» "the basic 
model includes the rotor inflow perturbation, since rotor- rotor interference 
has an important role in the dynamics. See references 9 and 10 for dlscusBlons 
of the modes represented by the eigenvalues given in tables 2 to 5* 

Generally it is concluded from the data in tables 2 to 5 that the 
quasistatic approximation is very good? that the inflow perturbation has 
a significant, and sometimes very important influence on the flight dynamics? 
and that the assumption of separate longitudinal and lateral motions is 
not usually valid, at least in terms of the eigenvalues, Rotor lag motion 
and other degrees of freedom may be required in the model? this is particularly 
true for the hingeless rotor, where the degrees of freedom have a large 
influence on the dynamics. 

CONCLUDING REMARKS 

A comprehensive aralysis of rotorcraft aeroelastlc behavior (reference l) 
has been applied to a number of classical problems in rotor dynamics. The 
results presented agree well with the solutions available in the literature. 

The satisfactory treatment of the classical problems such as these is a basic 
prerequisite for any analysis which is to be applied to the more complex 
problems of rotorcraft dynamics. In addition, it has been possible to 
examine the influence of a number of features of the analytical mode, including 
the effect of inflow perturbations on rotor flap response and helicopter 
flight dynamics, the characteristics of the ground resonance instability, and 
the role of the rotor degrees of freedom in helicopter flight dynamics. 
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1 , Rotor and support parameters for the ground 
resonance calculations of figure 9. 


Rotor I 

Radius 5.3 m 

Number of blades 4 

Hinge offset {% radius) 5*9 

lag frequency, (per rev) ,31 

Blade mass 46.3 kg 

Blade first moment of inertia II6 kg-m 

p 

Blade second moment of inertia 381 kg-m ' 


Support s 


mode 

natural 

frequency 

Hz 

modal 

mass 

kg 

modal 

damping 

N/m/sec 

long, balance 

1.76 

?,7000 

36500 

lat, balance 

E.16 

11000 

6900 

long, strut 

3*04 

24500 

38000 

lat, strut 

3.68 

18000 

36500 




- 1 - 



Table 2 . Small Articulated Rotor Helicopter Plight Bynzuaics 


I 

Co 

I 


a 

ts 

V- ♦ 

o 

o ^ 
o ^ 
3 « 
pa 


A . Hover 


1 . Basic model 

2 . Quasistatic rotor 
model 

3 . Including rotor lag 
motion 

4 . Including inflow 
perturbation 

5 . Coupled longitudinal 
and lateral dynamics 

B. V = 100 knots 


1 . Basic model 

2 . Qxiasistatic rotor 
model 

3 . Including rotor lag 
motion 

4 . Including inflow 
perturbation 

5 . Coupled longitudinal 
and lateral dynamics 


Longitudinal modes 


T 

^2 

ti 

2 

ti 

12.90 

3 . 9 ? 

1 .769- 

.903 

12.90 

4.00 

.77? 

.913 

12.90 

3.95 

.769 

.900 

12.90 

3.96 

1.170 

-B 93 

15 .^ 

1.957 

.700 

.816 


Longitudinal modes 


T 

ti 

! T ti 

i 2 

31 M 

15.65 

2.32 

.608 

31.42 

15.55 

2.34 

.619 

31.^1 

15.71 

2.32 

.610 

32.70 

i 

14.71 

2.37 

.645 

30.B6 

34.39 

2.26 

.610 



Lateral 

modes 


T 



ti 

p. 

11.32 

9.93 

.388 

1.137 

11.35 

10.00 

.399 

l.l 4 lf 

11.36 

10.14 

.386 

1.166 

11.12 

9.30 

.403 

1.238 

7.8it 

tl-15.34 

2 

T » 19.33 
ti = .600 

2 


Lateral modes 


T 

ti 

2 

T 


20.30 

2.06 

1.552 

.im 

20.31 

2.07 

1.558 

.484 

24.01 

1.92 

1.563 

.494 

20.3P 

2.03 

1.582 

.501 

20.21 

1.^1 

l. 5 fe 

.487 



Table 3- I»arge Articulated Rotor Kelit^pter Plight Dyaaaics 


A , Hover 


Longitudinal nodes Lateral nodes 



T 


ti 

2 

ti 

2 



% 

> 

H 

1 . Basic model 

19.93 


.535 

U9^ 

11.30 

10.?« 

.381 

2.23 

2. Onasistatic rotor 
model 

19.86 


.591 

1.599 

11.23 

IC.'^ 

.Tl 

2.2? 

\ 

3. Including rotor lag 
motion 

19.90 

5.3S 

.522 

1.519 

11.51 

10.72 

.379 

2.13 

Incliiding inflow 
perturbation 

19.92 

5.^1 

l,0?h 

1,597 

11.07 

18.32 

.i^36 

2.67 

5. Coupled longitudinal 
and lateral dynamics 

23.33 

^.17 

.56^?' 

1.372 

9 . 39 ^ 

tx^7.B3 

z 

.371 

2.35 


B. V - 100 knots 


Longitudinal modes lateral rodes 




T 

^2 

^ tl 

2 

tl 

2 

T 

tl 

2 

tl 

2 

4 

1 . 

Basic model 

33.80 

5.37 

.373 

2.72 

3.5? 

1.099 

.625 

5.59- 

2 . 

Quasistatic rotor 
model 

33.52 

5.40 

.408 

2.75 

3.60 

1.118 

.610 

5.57 

3. 

Including rotor lag 
motion 

3^.36 

5.37 

. 3?0 

2.70 

3.57 

1.143 

.586 

6.58 

9 ^. 

Including inflow 
perturbation 

39-.78 

5.23 

.471 

2 .88: 

3.B1 

i.326 

.607 

1 

5.58 

1 

5. 

Coupled longitudinal 
and lateral dynamics 

40.71 

5.^ 

T = 

tl = 
z 

19.32 
= .471 

3.63 

1.166 

1 

T = 48.43 

tl “ 3.38 
z 


*fe.l5le 4* Soft-Inplane Hingeless Ho tor Helicopter Flight IJynanics 


A • Hover 


B. 


Longituciinal aodes lateral Wides 



T 

*2 

ti 

- ■=’ i 

ti 

2. 

T 

_ ^ % 


*2 

1 . Basic model 

17.04 

4.21 

.639 

1.393 

5 .^ 

3.31 

.339 

1.2^ 

2 . Qi\asistatlc rotor 
model 

17.03 

4.23 

.648 

1.398 

6.03 

3.79 

.340 

1.^3 

3. Including rotor lag 
motion 

17.06 

4.20 

.63« 

1.336 

5.99 

3.60 

.330 

1.232 

4^. Including rotor lag 
and torsion motion 

17.0? 

3.65 

.783 

2.33 

6.33 

16.48 

.374 

.939 

5. InclTiding inflow 
perturbation 

17.95 

6.28 

I. 0 F 7 

1.0*^ ' 

7.12 

3.44 


1.366 

6. Coupled, longitudinal 
ar»i lateral dynamics 

23.56 

3.65 

.661 

1.726 

6.05 

2.84 

.354 

1.609 


V - 100 fcnots 

Longitudinal modes Lateral modes 




I 


T 

ti 

p 

T 

•> 

*> 

*2 

1. 

Basic mod'll 

63.19 

19.16 

3.38 

.77? 

1.859 

.785 

.611 

7*62 

2. 

Qu^istatic rotor 
model 

63.23 

19.09 

3.40 

.79> 

1.856 

.781 

.630 

?-65 

3 . 

Including rotor lag 
motion 

63.49 

19-28 

3.37 

.776 

1.861 

.792 

.602 

7.50 

4. 

Including rotor lag 
and tension motion 

36.30 

15.34 

3.02 

1.1^ 

1.785 

,?2^ 

» . 

21.65 

= 16.04 

5 . 

Including inflow 
perturbatJ >n 

70.92 

18.60 

3.33 

.848 

i .910 

.822 

.569 

8.34 

6. 

Coupled longitudinal 
and lateral dynamics 

254.3 

tg= 13.89 

3.42 

.863 

1.852 

.??3 

.540 

2 
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TaT)le 5 » Tandem Ho tor Helicopter Flight Dynamics 


A . Hover 


1 , Basic model 

Z* Quasistatic rotor 
model 

3 . Including rotor lag 
motion 

4 , Without inflow 
perturbation 

5 • Coupled longitudinal 
and lateral dynamics 


B . V =» 50 knots 


1 . Basic model 

2 . Quasistatic rotor 
model 

3 . Including rotor lag 
motion 

4 . Without inflow 
perturbation 

5 * Coupled longitudinal 
and lateral dynamics 


Longitudinal modes 


T 

^^2 

ti 


25.04 

27.37 

.675 

1.155 

25.06 

27.62 

.691 

1.172 

24.99 

27.10 

.671 

1.152 

26.15 

51.25 

.463 

.720 

25.06 

11.97 

.640 

1.155 


Longitudinal modes 


T 

*2 

tjL 

. . .51 

tl 

26.40 

4.00 

,322 

3.15 

26.31 

4.00 

.341 

3.16 

26.1? 

4.00 

.331 

3.13 

27.43 

98.'“'' 

. 44-7 

.711 

29.13 

3.83 . 

.322 

3.19 


lateral modes 


T 

*2 

tl j 

2 1 

tl 

2 

12.78 

4.75 

.759 ! 

15.01 

12.74 

4.82 

.766 

15.02 

12. 78 

4.70 

.750 

14.89 

12.78 

4.75 

.759 

14.37 

12.78 

4.72 

.776 

8.75 


Lateral modes 


T 

*2 

% 

2 

% 

10.49 

6.86 

.713 

9.17 

10.54 

6.84 

.716 

9.20 

10 .^ 

5 . 3 « 

.768 

10.46 

10.47 

6.83 

.714 

9.31 

10.62 

5.23 

.711 

4.23 



Table 5 


Concluded 


C. V = 100 knots 


1. Basic raodel 

2. Quasistatic rotor 
model 

3. Including rotor lag 
motion 

h, Without inflow 
pert’jrhation 

5. Coupled longitudinal 
and lateral dynamics 


Longitudinal modes 


T 

*2 

ti 

2 

34.05 

2.85 

.26? 

3.04 

33.94 

2.84 

.284 

3.05 

33.75 

2.82 

.271 

3.01 

25.48 

21.73 

.346 

1.071 

12.91 

2.10 

.266 

3.35 


Lateral modes 


T 

to ! 

j 

ti 

2 

ti 

2 

15.07 

7.11 1 

1.295 

1.295 

15.15 

6.36 

1.130 

1.422 

11.12 

3.88 

.937 

4.97 

15.05 

7.15 

1.298 

1 . 29 ^ 

48.62 

4.36 

T « 15.94 

ti 1.001 

2 



No infiow dynomics 
Qua$lstatic inflow 
With inflow dynamics 


Figure 1. Magnitude and phase of coning response to collective inputs 
( ) for a hovering, three-bladed articulated rotor 
(flap frequency >? = l.O), showing the influence of the 
inflow dynamics. 
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io/Sl 


Figure 2, Magnitude and phase of lateral tip-path-plane tilt 
response to lateral control plane inputs 
for a hovering, three-bladed hingeless rotor (flap 
frequency ^ = 1.15)f showing the Influence of the 
inflow dynamics; 
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Figure 3* Magnitude and phase of longituflinal tip- path-plane 

tilt response to lateral control plane inputs ((3ic/dic.) 
for a hovering, three-hladed hingeless rotor (flap 
frequency = 1 . 15 )» showing the influence of the 
inflow dynamics. 
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Figure 4. Root locus of forward-flight flapping stability for a 

single articulated rotor blade ( V = 10, flap frequency 
= 1.0, no pitch/flap coupling). 
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xa/C =0 

Xa/C = .05 


.10 


05 


of an articulated 





Figvire 6, Pitch/flap flutter stability toundaries of a hovering 


articulated rotor, with and without the first flatwise 
"bending mode. 











Figure 8. Comparison of blade lag damping at the ground resonance 

stability boundary, as calculated by the Deutsch criterion 
and using the dynamics analysis of reference 1. 
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0 100 200 300 400 

rpm 


i'^igure 9. Ground renoiiance stability calculations for a four-bladed 
articulated rotor on a -winr! tunnel strut and balance 
system! lamping ratio of the support modes as a function 
of rotor speed, for three levels of blade lag damping. 



